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Abstract. We analyzed ionized gas motion and disk orientation parameters for 15 spiral galaxies. Their velocity 
■ fields were measured with the Ha emission line by using the Fabry-Perot interferometer at the 6m telescope of 

SAO RAS. Special attention is paid to the problem of estimating the position angle of the major axis (PAq) and 
the inclination (i) of a disk, which strongly affect the derived circular rotation velocity. We discuss and compare 
' different methods of obtaining these parameters from kinematic and photometric observations, taking into account 

the presence of regular velocity (brightness) perturbations caused by spiral density waves. It is shown that the 
commonly used method of tilted rings may lead to systematic errors in the estimation of orientation parameters 
(and hence of circular velocity) being applied to galaxies with an ordered spiral structure. Instead we recommend 
using an assumption of constancy of i and PAo along a radius, to estimate these parameters. For each galaxy 
of our sample we present monochromatic Ha- and continuum maps, velocity fields of ionized gas, and the mean 
Oh' rotation curves in the frame of a model of pure circular gas motion. Significant deviations from circular motion 

with amplitudes of several tens of kms -1 (or higher) are found in almost all galaxies. The character and possible 
5_j ■ nature of the non-circular motion are briefly discussed. 
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, 1. Introduction circular motion may be different. Small-scale gas motion 

is usually related to star formation sites (HII complexes), 
Most of the available data on the kinematics of ionized gas whereag large _ scale velocity per turbations are caused by 
in disk galaxies is obtained with long-slit spectrographs, waye spir£j pattem &nd/or central oyal (bar) impacts Qn 
allowing us to construct mean rotation curves under an ^ e g aseous disk 
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assumption of pure circular rotation. But the circular ro 
tation model, being useful for determinations of the global 
properties of a galaxy, is only a zeroth order approxima- 
tion to the real and rather complex dynamical picture. 



At the end of the 80s, a series of observations with 
a long-slit spectrograph was undertaken at the 6m tele- 
scope of the SAO RAS, which have revealed significant 

. non-circular motion with amplitudes of 30 — 100 kms -1 

Fven if there is no active nucleus or tidal perturba- . . , , . /Ar . , , noo mnr>\ 

f . in ordinary spiral galaxies (Atanasiev et al. 1988 - 1992). 

tion from a nearby companion to affect gas kinematics t, , , , , . r ,. , . 

But tne most complete information on tne gas kmemat- 
m a galaxy, the observed non-circular motion may have . > , , n , n , • r n r,T-> 

,7 , i .... ics m galactic disks can be acquired only by using full 2D 

amplitudes up to 30 — 50 kms (or even higher in some , n , , , • , , , , . , . . , , . n 

^ , v ° velocity fields which may be obtained m the optical spec- 

cases), being therefore supersonic. The nature of non- , , , r . , n -r\\ ± 

n ° ^ tral range by methods of panoramic (2D) spectroscopy. 

* Based on observations collected with the 6m telescope of The method based 011 the Fabry-Perot interferometric ob- 

the Special Astrophysical Observatory (SAO) of the Russian servations was used by several groups of observers to get 

Academy of Sciences (RAS), operated under the financial sup- the velocity fields for large samples of different morpho- 

port of the Science Department of Russia (registration number logical types of galaxies (see f.e. Schommer et al, 1993, 

01-43). Amram et al. 1995). The largest project of this type is 
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the GHASP survey, carried out at the Haute-Provence 
Observatory, which should provide a homogeneous sam- 
ple of 2D velocity fields of about 200 spiral and irregular 
galaxies. This survey aimes to compare galaxies in dif- 
ferent environments, different stages of evolution and at 
different redshifts (Garrido et al. 2002, 2003). Unlike the 
works mentioned above, the main goal of our project is 
to study the regular deviations of gas motion from circu- 
lar rotation in spiral galaxies, which may be caused by 
spiral density waves - in parallel with the routine task 
of obtaining the rotation curves and mass distributions. 
For this purpose we use the Scanning Fabry-Perot inter- 
ferometer (IFP) of the 6m telescope. We have obtained 
high-resolution line-of-sight velocity fields of the ionized 
gas over a large field of view (about 3') for 42 spiral galax- 
ies, 1 both barred and non-barred, mostly of Sb-Sc types, 
by measuring the Ha (and in some cases also [Nil]) emis- 
sion line. This sample may be by no means considered 
as a complete one to some limiting magnitude or angular 
size. It just represents spiral galaxies of different types 
with the angular diameter of 2' < D25 < 6' (for the 
most of them) having various types of spiral patterns - 
from certain grand-design to flocculent ones. The results 
of the velocity field analysis for some of the galaxies of 
our sample have already been published (Fridman et al. 
1997, 1998, 1999, 2001a,b,c; Fridman & Khoruzhii, 2003). 
In this work we present the results of observations and 
the rotation curves for 15 galaxies of the sample. In paral- 
lel with the Doppler velocities measurements, brightness 
maps in Ha and near-Ha continuum were also obtained 
for each galaxy. 

In principle, the Doppler velocity field contains infor- 
mation both on the rotation and on non-circular gas mo- 
tion, including the ordered wave motion related to a spiral 
density wave (if it exists in a chosen galaxy). The prob- 
lem to be solved is the development of reliable methods 
which would enable us to extract and analyze these mo- 
tions separately. Historically, the first attempts at such an 
approach were based on drawing a "mean" rotation curve, 
and on the consequent extraction of non-circular motion as 
residuals left after the subtraction of the circular rotation 
component from the observed velocity field (it was done 
for the first time by Warner et al. (1973) for M 33 and 
Pence (1981) for NGC 253). The attempts to relate the 
non-circular velocities found in this way to the observed 
spiral structure were unsuccessful, and this was treated 
by the authors as a failure to find non-circular motion of 
wave nature in these galaxies. 

On the other hand, it was known from earlier papers 
that there exists a correlation between the velocity pertur- 
bations and the pattern of spiral structure in some galaxies 
(see for example Rubin et al. (1980) for optical rotation 
curves of galaxies or Rots & Shane (1975) for a HI velocity 
field of M 81). If that is the case, can the negative result 
obtained for M33 and NGC 253 be accounted for by some 

For the full list, see 
http:/ /www. sao.ru/hq/lsfvo/vortex/vortex. html 



peculiarities in these spiral galaxies? Indeed, although 
M 33 has two main spiral arms, it is really a multi-armed 
galaxy (Sandage & Humphreus 1980). In turn, NGC 253 is 
a starburst galaxy with a rather high disk inclination, in- 
appropriate for tracing spiral arms. However, later works 
by Sakhibov and Smirnov (1987, 1989, 1990) made it ev- 
ident that the main problem of these early approaches 
is that independent determination of the rotation curve 
and of the residual velocities related to a density wave 
is not possible even for an "ideal" two-armed grand de- 
sign galaxy. This conclusion was based on the analysis of 
the line-of-sight velocity fields of some galaxies by apply- 
ing a Fourier-series expansion of velocity along the galac- 
tocentric azimuthal angle. The authors mentioned above 
showed that the even component of the first Fourier har- 
monic includes contributions from both the rotational ve- 
locity and the motion related to the density wave in the 
disk plane. Such interference makes it necessary to work 
out a special method to separate the circular and non- 
circular motion consistently. This method was proposed 
and described in detail in the papers by Lyakhovich et al. 
(1997), Fridman et al. (1997, 2001). These authors also 
showed that, in general, ignoring the influence of system- 
atic (density-wave related) non-circular motion onto the 
line-of-sight velocities may prompt systematic errors in es- 
timations of orientation parameters and lead to a wrong 
conclusion about their radial variations, even in the case 
of their constancy. 

A similar problem appears when the photometric data 
is used to determine i or PAq of a disk. In this case the 
isophote distortion due to spiral arms may also mimic the 
radial variations of orientation parameters. To avoid this 
difficulty the outer disk isophotes are traditionally used 
where the contrast of spirals is usually low and may be 
neglected. But in some cases these regions may not be 
covered by the field of view, or are located too far from 
the galactic center so that their orientation differs from 
that of the inner disk (say, due to a disk warp). In this 
situation the correct estimation of orientation parame- 
ters needs a complex self-consistent analysis of brightness 
and/or of line-of-sight velocity distributions, taking into 
account both non-circular motion and brightness asym- 
metries, but this way is rather cumbersome for practical 
use. A simplified method should be chosen which may give 
an acceptable accuracy of orientation parameters. As it 
will be shown below, this problem is not trivial, and the 
use of different methods may lead to diverging results. 
Note that the accuracy of estimation of inclination and 
position angle of a galaxy influences not only the rotation 
curve obtained from observations. It may be a key factor 
which determines the possibility and the reliability of the 
evaluation of residual velocity field and of some refined 
parameters which, for example, may be used to distin- 
guish chaotic and regular trajectories of gaseous clouds 
(Fridman, Khoruzhii & Polyachenko 2002). 

Section 2 contains descriptions of our observations and 
data reduction procedures. In Section 3 we discuss various 
approaches to the determination of orientation parameters 
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from the kinematic and photometric data and compare 
two independent methods, the combined use of which en- 
ables us to solve the problem with good accuracy. The 
parameters we found for 15 galaxies are used to obtain 
their velocity curves in a model of pure circular motion 
of gas. The discussion of individual galaxies is given in 
Section 4. Different types of non-circular gas motions are 
discussed in Section 5. The main results of this work are 
summarized in Section 6. 

The detailed analysis of the systematic motions re- 
lated to the density waves, based on the method described 
above, will be given in subsequent papers. 

2. Observations and data reduction 

The velocity fields of the chosen galaxies obtained in the 
Ha line cover significant parts of their discs, which enables 
us to use them to compare the estimates of the orien- 
tation parameters from both kinematic and photometric 
data and to obtain the rotation curves of the galaxies. 

The observations were made in 1995-2000 with the 
Scanning Fabry-Perot Interferometer (IFP) attached to 
the reducer at the prime focus of the 6m telescope of the 
SAO RAS; the equivalent focus distance of the system is 
F/2.4. The description of the interferometer can be found 
in Dodonov et al. (1995). 

We limit the required spectral range by using a narrow 
interference filter with the FWHM = 10 - 20 A and a 
transmittance peak near a redshifted Ha line of a galaxy. 
When observing some galaxies (IC 1525, NGC 1084, etc.), 
the emission lines [Nil] A 6548, 6583 from the neighbouring 
interference orders weakened by the wings of the filters 
also get into the passbands. Fortunately, when using any of 
our interferometers, the interference orders were overlaid 
in such a way that the gap between Ha and [Nil] was not 
less than half of the full spectral range, so the presence of 
[Nil] in the spectra did not prevent exact measurement of 
Ha, excluding the case of the central part of NGC 5371 
which is a LINER possessing rather broad emission lines. 

During various observational runs we used two IFP by 
the Queensgate company which are operating in 235th and 
501st interference orders (at a wavelength of 6562.78 A) 
with the finesse of 8-12. The parameters of the interferom- 
eters are given in Tabled The number of spectral chan- 
nels (or steps of scanning) for various objects was 24 or 
32. Also, various types of panoramic detectors were used 
for different observational runs. In 1995-1996 we used a 
2D photon counter (IPCS) with a format of 512 x 512 
pixels, in 1997 - a CCD with a format of 1040 x 1160 
pixels, from 1998 to the present - the CCD TK1024 with 
a format of 1024 x 1024 pixels. In all cases the observa- 
tions were made with binning of 2 x 2 pixels to increase 
the S/N ratio and to economize the CCD readout time. 
The data are presented as "data cubes" consisting of 24 
or 32 object images obtained with varying gaps between 
the interferometer plates. 

For the different types of detectors we used different 
regimes of scanning. When observing with the IPCS, we 



Table 1. IFP parameters. 



Fabry-Perot interferometers 
FP235 FP501 



An order (at Ha line) 


235 


501 


Interfringe (AA) 


28A(1270kms- 1 ) 


13A(600 kms" 1 ) 


Resolution (S A) 


3A(130 kms" 1 ) 


1.2A(130 kms" 1 ) 


Spectral sampling: 






(24-channel mode) 


1.20A(53 kms" 1 ) 


0.55A(25 kms" 1 ) 


(32-channel mode) 


0.90A(40 kms" 1 ) 


0.40A(20 kms" 1 ) 



accumulated an image in every channel during 15-20 sec, 
and then the channel was switched forward. After pass- 
ing through all the channels, the scanning cycle was re- 
peated. This approach allowed us to overcome variations 
of atmospheric transparency and of seeing quality during 
total exposure. When observing with a CCD, each channel 
was exposed only once, but to properly take into account 
variations of atmospheric transparency, the odd channels 
(1,3,5...) are exposed first, and then the even channels 
(2,4,6,...). 

To correct spectra for the phase shift, we observe an 
emission line of the calibrating lamp. When observing with 
the IPCS, the calibration cubes were exposed before and 
after each object exposure; when observing with a CCD 
- at the beginning and at the end of the night. To cor- 
rect the data for the modulation introduced by the order- 
separating filter and for the flat field, the telescope main 
mirror cup or the wall of the dome lit by a continuum light 
source are observed in every channel. 

Reduction of the data obtained with the IPCS was 
performed in a standard way with the software ADHOC 
(Boulesteix 1993). It included a construction of the phase 
map by using the cube of the comparison spectrum, 
flat-fielding, phase-correcting the object cube (wavelength 
scale calibration), and sky spectrum subtraction based on 
analysis of detector areas free of the object emission and 
of its ghost images. To make an optimal filtration of the 
data and to increase the S /N ratio in the regions with weak 
emission, we smoothed each spectrum in the spectral do- 
main by a ID Gaussian with FWHM = 2-3 channels and 
in the space domain by a 2D Gaussian with FWHM close 
to the seeing estimated from the measurements of nearby 
stars. 

The primary reduction of the data obtained with 
CCDs is performed with the software developed by us 
and described in Moiscev ( 2002a 2002b ). It includes bias 
subtraction, flat-fielding, cosmic hit cleaning, and mu- 
tual shifts of frames to overlay reference stars. The same 
stars are used to correct the datacube, if necessary, for 
the variations of atmospheric transparency and of see- 
ing quality. Positions of the comparison spectrum line 
are determined by fitting its profile by a Gaussian or 
by a Lorenz profile; this modification has allowed us 
to achieve somewhat higher accuracy in the phase map 
construction with respect to the previous results ob- 
tained with ADHOC. After night-sky spectrum subtrac- 
tion and correction for the phase shift, further reduc- 
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tion is performed by using the software ADHOC and a 
new version of this software for WINDOWS, ADHOCw 
I http:/ /www-obs. cnrs-mrs.fr/ADHOC/adhoc. html). 

The journal of observations is given in Table |2 where 
for each galaxy we list its morphological type according 
to RC3, the date of the observations, the detector and the 
interferometer types which were used in each particular 
observation, pixel size after binning, total exposure time 
in each channel, seeing quality, and spatial resolution after 
all smoothing procedures estimated by using completely 
reduced star images in the cubes. 

For each galaxy we have constructed a velocity field 
of the ionized gas by measuring barycenters of the Ha 
emission line in each spectrum (for one-component line 
profiles), as well as monochromatic images in Ha and red 
continuum images in the narrow passband near Ha. If the 
emission-line profile is symmetric, the accuracy of the ve- 
locity measurements is restricted only by the accuracy of 
phase map constructing and is about 10 — 15kms _1 for 
the observations with the IPCS and about 5 kms -1 for the 
CCD data. However in some cases (see below) the emis- 
sion lines are clearly double-humped or possess strongly 
asymmetric profiles which is evidence for the presence of 
several subsystems of ionized gas clouds on the line of 
sight. In these cases we have performed a multi-component 
Gauss analysis of such profiles to separate the kinematic 
components. The resulting errors in the "main compo- 
nent" velocity measurements are about 15 — 20 kms -1 for 
FP235 and 5 — 10 kms -1 for FP501 observations, which is 
similar to the errors of barycenter-based velocities in the 
case of a symmetrical line profile. To construct a veloc- 
ity field in the regions with multicomponent profiles we 
use only those velocities which correspond to the "main' 
component" . In practically all the cases its value appears 
to be nearest to the velocities found in the neighbouring 
regions with a one-component shape of an emission-line 
profile. With the exception of a complex situation in the 
central region of NGC 1084, such an approach has not cre- 
ated any velocity discontinuity. The examples of the emis- 
sion line spectra in the regions with the "abnormal" pro- 
files are plotted in Fig. ^ Double-peaked profiles may be 
found in the circumnuclear regions (NGC 23, NGC 1084, 
NGC 1134, NGC 2964, NGC 6643), as well as in the outer 
disks (IC 1525, NGC 1084, NGC 3893). 

The most complex case is represented by the circum- 
nuclear regions of galaxies with broad stellar Ha absorp- 
tion lines in the spectra. Unfortunately, our spectral range 
is too narrow, and the IFP order overlapping often pre- 
vents a proper correction for this absorption feature. We 
met such situations in some regions in the centers of 
IC 1525, NGC 4136, NGC 4414, and NGC 5371. We have 
tried to use a double-Gaussian (emission + absorption) 
model to calculate the velocities in these regions. However 
the errors of the velocity measurements were too large 
(~ 50 — 70kms~ 1 ). Therefore we have excluded these re- 
gions from the resulting velocity fields (see the descrip- 
tions of the individual galaxies in Section 0}. 



3. Results 

3.1. "Kinematic" orientation parameters of galactic 
disks 

After gas velocity fields of galactic disks became available 
(first from HI radio line observations, later from optical 
line measurements), several ways were proposed to solve 
the problem of restoring rotation curves of galaxies from 
the obtained data. In most cases the basic assumption is to 
consider deviations of the observed velocities from those 
expected in a pure circular rotation picture as random 
ones, which have a "noisy" character: 

yobs _ ymod _|_ q, ^ 

V mod (r, <p) ee V sys + U rot (r)cos<^sinz, (1) 

where ^ is a random deviation of the observed velocity 
from that predicted in the model; V sys is the systemic line- 
of-sight velocity of the galaxy; V ro t is the model circular 
velocity of the gas in the plane of the galaxy which de- 
pends on the galactocentric radius r; tp is the galactocen- 
tric azimuthal angle measured with respect to the major 
axis; i is the inclination (the angle between the galactic 
plane and the sky plane). 

But even this very simplistic approach leaves the prob- 
lem of determining the rotation curve to be rather compli- 
cated. In general, there are five parameters which should 
be found or taken from elsewhere in order to extract 
the rotation velocities from the line-of-sight velocity field. 
They are: the systemic velocity of a galaxy V sys , the posi- 
tion of the center of mass xo^yo, the position angle of the 
line of nodes PAq and the disk inclination angle to the 
line of sight i. 

In early papers (Warner et al. 1973, Pence 1981) these 
parameters were assumed to be constant over the entire 
disk. In this case they could be found either from the 
analysis of photometric data (with the exception of V^ ys ) 
(Pence 1981), or by minimizing a sum of squared differ- 
ences between the measured and model-predicted veloci- 
ties for a pure rotation model (Warner et al. 1973, Pence 
1981). In this case the rotation velocity at a given galac- 
tocentric radius was calculated as a weighted mean of all 
points lying in the plane of the disk within a ring of ra- 
dius r and of width dr. Even in these early works it was 
revealed that different methods give discrepant results, 
which was ascribed to the influence of non-circular mo- 
tion (Pence 1981). Note that a similar effect could also 
be caused by non-planar motion. As it was mentioned by 
van Moorsel and Wells (1985), the Warner et al (1973) 
method is of a trial-and-error nature, which makes it diffi- 
cult to determine (cite) "what errors are to be assigned to 
the results" . These authors proposed a modification of the 
ring method assuming a certain functional approximation 
of the rotation curve under the best-fitting orientation pa- 
rameters; after that they solved the system of many equa- 
tions for the total number of unknowns by a least-squares 
algorithm (without using any weighting based on the posi- 
tion inside of a galaxy) . Unlike van Moorsel & Wells, below 
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Table 2. Log of IFP-observations. 



Name 




Type(RC3) 


Date 


Detector 


IFP 


px size 


T exp , s 


Seeing 


Beam 


IC 


1525 


SBb 


04 Dec 99 


CCD 


FP235 


0'.'68 


24 x 180 


2'.' 5 


3'.' 2 


NGC 


23 


SB(s)a 


29 Aug 98 


CCD 


FP235 


0'.'68 


24 x 180 


2'.'0 


2'.' 8 


NGC 


615 


SA(rs)b 


29 Aug 98 


CCD 


FP235 


0'.'68 


24 x 180 


2'.'5 


3'.'8 


NGC 


972 


SA(rs)b 


21 Feb 98 


CCD 


FP235 


0'.'68 


24 x 180 


2'.'6 


3'.'5 


NGC 


1084 


SA(s)c 


26 Oct 95 


IPCS 


FP501 


Of! 92 


32 x 225 


2'.'0 


3'.' 5 


NGC 


1134 


S..? 


23 Feb 98 


CCD 


FP235 


0'.'68 


24 x 150 


2'.'4 


2'.'9 


NGC 


2964 


SAB(r)bc 


01 Dec 99 


CCD 


FP235 


0'.'68 


24 x 180 


2'.'0 


2'.'5 


NGC 


3583 


SB(s)b 


28 Feb 00 


CCD 


FP235 


0'.'68 


32 x 180 


2'.'8 


3'.'5 


NGC 


3893 


SAB(rs)c 


10 Dec 96 


CCD 


FP501 


0'.'90 


24 x 180 


2"5 


3'.'1 


NGC 


4100 


R'SA(rs)bc 


21 Mar 96 


IPCS 


FP235 


0'.'77 


24 x 220 


2'.'5 


3'.'5 


NGC 


4136 


SAB(r)c 


21 Feb 98 


CCD 


FP235 


l'/36 


24 x 180 


2'.'5 


4'.'0 


NGC 


4414 


SA(rs)c 


21 Feb 98 


CCD 


FP235 


0'.'68 


24 x 120 


2'.'0 


3'.'0 


NGC 


4814 


SA(s)b 


23 Feb 98 


CCD 


FP235 


Of! 68 


24 x 120 


2'.'4 


3'.'1 


NGC 


5371 


SAB(rs)bc 


23 Mar 96 


IPCS 


FP235 


0"77 


24 x 460 


4'.'0 


5'.'5 


NGC 


6643 


SA(rs)c 


16 May 97 


CCD 


FP235 


0'.'45 


24 x 360 


l'/5 


2'.'0 



we do not use any parameterization of the rotation curve, 
applying the least-square algorithm to the model of a thin 
circularly rotating disk where the velocity of rotation of a 
given ring is considered as the unknown. 

The parameters PAq and i estimated in the early 
works were found in some cases to be inapplicable to the 
outermost parts of galaxies, where a distortion of isoveloc- 
ities is often seen (see for example Warner et al. 1973 for 
M 33, or Rogstad et al. 1979 for NGC 300). To improve 
the procedure without rejecting a circular motion approx- 
imation, it was proposed (Begeman 1989) to consider the 
values of PAq and i (or in general, all the parameters in- 
volved) as functions of r. This approach is called "a tilted- 
rings model" . In this model it is assumed that the disk of 
a galaxy may be described by a set of concentric rings 
with fixed values of V sysi Xo, z/q, i, and PA within each 
of them. In other words, each ring is allowed to have ar- 
bitrary circular velocity and orientation parameters. The 
values of V sys , xq, and j/o are often taken as the same for 
all the rings. This approach is widely used to determine 
rotation curves. 2 

Variations of the parameters are allowed not only at 
the periphery of a disk, that may be physically advocated 
because of a possible warp, but at all radii. It is worth 
mentioning that although the warping of gaseous (and 
sometimes stellar) disks of galaxies is a well-known fact, 
significant variations of inclination usually take place in 
the outermost parts of galaxies and, as a rule, beyond 
the regions covered by optical spectral measurements. It 
seems that the parameter variations obtained by formal 
use of the tilted- rings model is often a consequence of non- 
circular motions which are not taken into account by this 
method. If the variations of either i or PAq (or xo, yo) 

2 Note however that in some papers (see Schommer et al. 
1993, Amram et al. 1995, Garrido et al. 2002, 2003) the authors 
use constant values of orientation angles PAq and i for each 
galaxy as the best fit parameters, which may be found in the 
frame of circular rotation model by the least-square method or 
by iterations. 



bear a systematic character, it may be an indication that 
the presence of a spiral density wave (or bar) should be 
taken into account (Lyakhovich et al. 1997, Fridman et al. 
1997). 

As an illustration, Figs.6|21a, b show the radial varia- 
tions of i and PAq, formally determined by the tilted- rings 
method for NGC 3893 - one of the galaxies of our sample 
with rather moderate inclination. The systemic velocity 
and center position were assumed to be invariable. 

In Figs. |21 c, d analogous dependencies are shown for 
an artificial (model) velocity field. This field was created 
as a sum of pure rotation, non-circular regular velocities 
expected in the simplified model of two-armed density 
waves, and a "noisy" non-correlated component imitat- 
ing the errors of the velocity measurements and the oc- 
casional small-scale velocity dispersion. For definiteness, 
the amplitude of the rotational component was taken to 
be 210 km s -1 , and of the regular non-circular velocities - 
10 — 20 kms -1 , and the root mean square deviation of the 
"noisy" component was taken as lOkrns -1 . Circles mark 
estimates of the parameters when regular non-circular ve- 
locities induced by the density wave are excluded. One can 
see that dependencies obtained for the model which in- 
cludes the density waves reproduces the radial variations 
of parameters formally found for the real galaxy rather 
well. It is also clear that when the density wave is elimi- 
nated from the model velocity field, the random motions 
may lead only to small parameter variations (by 1-2 de- 
grees). It seems that the large variations of parameters 
in the tilted-rings model (Figs. [3 a, b) may be naturally 
explained by the presence of a density wave. 

The natural starting point for determining the orien- 
tation parameters is an assumption that they remain con- 
stant over all the optical part of a disk. But this in itself 
is not enough to find the angles we seek in a correct way if 
regular velocity perturbations in density waves take place. 
The main problem is that the latter introduce systematic, 
non-occasional deviations from the pure circular motion, 
that do not allow the use of the representation of the mea- 
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Fig. 1. The examples of the IFP-spectra with complex structure of the Ha emission- line profiles (see Section 0] for 
the detailed description of the individual galaxies). The thick line shows the spectrum extracted for a region with the 
cylindrical coordinates (r, PA) given at the top-right corner of the each frame. These profiles reveal different kinds 
of peculiarities, such as blue/red wings (see c, d, e, f, h, i and 1), double-humped structure (see g, j), f/a-absorption 
distortion (see b and k), and the [Nil] emission line interference (see a, g, and h). The thin line reproduces the spectrum 
of the nearest galactic region with the "normal" (symmetrical) shape of the emission line profile. The arrows mark 
the position of the [Nil] (A6548A or A6583A) line, intruding fr om the neighboring interference order. Note that the 
wavelength scale corresponds only to the Ha line velocities. 



sured line-of-sight velocity in the form JJJ. Due to the 
same reason, the orientation parameters calculated for a 
given ring may be shifted in a systematic, non-occasional 
way from their real values. It also renders it useless to seek 
the most probable values by averaging the data obtained 
for different rings, as it takes place often (see for example 
Carignan et al. 1988, Begeman 1989, Beauvais & Bothun 
1999, 2000, 2001). 

In order to not resort to the complex and laborious 
method of constructing a self-consistent model of non- 
circular motion in the disk, which being applied to an 
observational velocity field may give the best fit param- 
eters (see for example Fridman et al. 2001a, 2001b), one 



may use a simplified approach. The idea was proposed by 
Fridman et al. (2001a), where a radial periodicity of mo- 
tion in the density wave leading to a periodical character 
of the wave contribution into the Doppler velocity compo- 
nent is suggested. Hence, if one is to seek the best values 
of the orientation parameters, one can expect that among 
the methods ignoring the density wave component of non- 
circular motion, the best method uses the velocity field 
as a whole, without estimating individual parameters for 
separate rings. In this case the contributions of systematic 
deviations from circular rotation in different parts of the 
disk are self-compensated, at least partially. 
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Fig. 2. Radial variations of the inclination angle (left) and position angle of the line-of-nodes (right) found by the 
tilted-rings method applied to the observed velocity field of NGC 3893 (upper diagrams) and to the artificial galaxy 
(below), in which the velocity field is a sum of pure rotation, non-circular regular velocities expected in the simplified 
model of density waves, and "noisy" velocity dispersion of about 10 km/s. Thin lines in the lower diagrams mark 
the values of parameters obtained under the assumption of unique orientation parameters for the whole velocity field. 
Circles are the estimates made for the model in absence of density waves. Thick lines are real values of parameters 
accepted in the model. 



Actually this approach (without discussing the role 
of non-circular motion) was first used by Warner et al. 
(f 973). Our basic principles of the determination of rota- 
tion curves remain the same as those used by Warner et 
al. (f973) aside from the different approach to weighting 
the data. 

A consistent and reliable way to estimate different 
components of the Fourier spectrum of the line-of-sight 
velocity field (and particularly the rotation velocity) was 
proposed and validated by Lyakhovich et al. (1997), 
Fridman et al. (f997, 2001a), Burlak et al. (2000). ft is 
based on the application of a least-square method to the 
expansion of the azimuthal distribution of the observed 
gas velocities into the truncated Fourier series. The rota- 
tional velocity is allowed to vary from one ring to another. 
This is the only parameter which runs along the radius in 
our model. 

For a given ring of radius ri and width dr a sum of 
squared deviations of the measured velocities from those 
expected in the model of pure circular rotation is: 

N, 

X 2 (n) = ^2 ( V ?'( r J' Vj) - ^sys ~ Vrot(n) COSlfj silll) 2 .(2) 
3 

Here the summation is performed over all N\ points 
of the ring where the line-of-sight velocity was estimated, 
and the weighting coefficients are assumed to be equal. 3 
Analogous expressions may be easily written for any choice 
of the weighting function. For given values of the center 
position, systemic velocity, inclination, and position angle 

3 We ignore a weighting function in this work. The use of 
cos <fi as a weighting function is hardly justified here because 
the parameters we are seeking are determined not from com- 
bining individual velocity measurements, but rather from the 
functional dependence of the line-of-sight velocities on the az- 
imuthal angle, the expression for which is supposed to be 
known. In this case each point has a weight which does not 
depend on ip. See also Begeman (1989) who wrote (cite): ".. In 
principle there is no objection to the use of equal weights, ex- 
cept when the uncertainties in the determination of the radial 
velocities are known to be different. (In general these uncer- 
tainties are not known.) " However he adds: "The use of non- 
uniform [azimuthal] weights [is] preferable. ..radial velocities in 
positions close to the minor axis carry less information about 
the underlying circular velocity than at positions close to the 
major axis. Also by this reason one should give more weight to 
positions close to the major axis" . Indeed the latter statement 
is correct only when the circular rotation is restored from the 
individual local velocity estimates, not from the comparison 
with the model dependencies expected for a flat disc. 



of the line of nodes the minimum of x2(n) is reached under 
dx^(n) / dV mt {n) = 0, that gives: 

yrot ( n ) = J_ EjCOBy>j(^fa, tfj)- Vsys) 

sin i Ylj cos2 <Pj 

For complete filling of the ring by the measurements, this 
expression is reduced to a familiar one for the reverse 
Fourier transformation. 

Under condition @ the quadratic sum of deviations 
in a ring takes the form: 

^(n) = ^{V j (r j ,<p j )-V Bya f - 

3 

r i 2 

— - 2 L - (4) 

For the whole disk we have: 

X 2 ( x o, Ho, i, PA, V sys ) =J^X 2 (n)- (5) 

i 

Some parameters contained in the brackets on the left 
side of Eq.©, (say, the center position and/or systemic 
velocity) may be considered as knowns. The others may 
be found by minimizing the expression JSJ. In practice 
we estimated V S y S and the center position as the best fit 
symmetry center of the observed velocity field by compar- 
ing the velocities of points situated symmetrically with 
respect to the center, and kept these values fixed. 

If the rotation curve has a steep gradient in the cen- 
tral part of a galaxy, it may be more convenient to use a 
slightly different version of the expression for a quadratic 
sum of deviations. It may be obtained in a similar way as 
written above, but starting from a constant angular veloc- 
ity for a given ring instead of a constant linear velocity. 
The resulting expression differs from by the replace- 
ment of all values of cos <pj onto Xj — rj cos (fj . 

For the velocity field of the model galaxy which in- 
cludes density wave perturbations, the parameters i and 
PAq found by the method described above are shown by 
thin lines in FigsElc, d. They coincide with the parameters 
of the model (thick lines) within a range of under one de- 
gree. It enables us to conclude that the method gives cor- 
rect estimates of the orientation parameters, even in the 
presence of systematic non-circular motions in a galaxy. 

"Kinematic" orientation parameters of our sample 
galaxies, obtained by this method, are presented in the 
third column of Table Their calculation was made by 
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minimizing the expression JSJ with fixed positions and sys- 
temic velocities of the center. Confidence intervals given 
in Table |3] were estimated from the formally determined 
errors of the mean values in a model of circular rotation 
and from the shape and width of the profile of x2(i, PAq) 
near its minimum. The obtained results are insensitive to 
the chosen width of rings until it is smaller than the ra- 
dial scale of velocity variations in a galaxy. The PAq is 
traditionally measured counter-clockwise from the North 
direction to the receding part of the major kinematic axis 
of the galaxy. 

Note that our method of interferometer data calibra- 
tion may introduce some small systematic errors in sys- 
temic velocity estimation (however these errors do not 
influence the measurements of relative velocity varia- 
tions in the line-of-sight velocity field). A comparison of 
our V S ys (the second column in Table OJ with the NED 
database values shows that the velocity difference is about 
10 — 50kms _1 for different objects - that is of the same 
order as the scatter of estimates obtained by different au- 
thors. 

The fourth column of Table|31contains the kinematially 
determined orientation parameters obtained by a formal 
averaging of values found for different rings (a classical 
tiltcd-rings model). Because of the non-random variations 
of the orientation angles of individual rings, it is incorrect 
to characterize the confidence intervals of these estimates 
by the formal value of the error in mean. For this reason, 
a root mean square deviation from the mean value was 
taken as a confidence interval. As expected, this method 
of averaging gives estimates that are not as reliable as 
those in the method we have preferred above. 

3.2. "Photometric" orientation parameters of galactic 
disks 

Photometric methods are widely used because they are 
much simpler and more understandable and do not de- 
mand analysis of the detailed velocity fields. The methods 
are usually based on estimating ellipticity and isophote 
orientation by fitting ellipses to the outermost isophotes 
which can be reliably traced (they are assumed to be 
undisturbed by disk warping or by internal absorption). 
Another, more sophisticated way of getting the results is 
to use azimuthal Fourier analysis of brightness distribu- 
tion of a galaxy in the plane of sky (Grosbol 1985) or in 
the plane of the galactic disk (lye et al. 1982). 

As in the case of kinematic methods, the main problem 
of the photometric approach is to take into account the 
deviations from the axial symmetry caused by spiral arms 
or a bar. It may be especially severe when a galaxy does 
not possess a radial zone wide enough where the bright- 
ness distribution stays axisymmetric. Ignoring this prob- 
lem may lead to large systematic errors in evaluating the 
orientation parameters. 

In a formal sense, the second Fourier harmonic may ap- 
pear as a result of erroneous orientation parameters which 



are used for the deprojection of an axially symmetrical 
disk. This harmonic also appears if the parameters used 
are correct, but non-symmetrical details such as spirals 
or a bar are present, so it is impossible to separate these 
factors within the frame of Fourier analysis. 

The only way to solve the problem is to use differ- 
ences in the phase and amplitude behavior of the second 
Fourier harmonic when different factors are responsible for 
its appearance. Particularly, the amplitude of the second 
harmonic caused by an error in the estimation of the in- 
clination grows as this error increases. If the inclination 
is underestimated, the line of maxima of the second har- 
monic would align with the major axis because of this 
error, and in the opposite case it would coincide with the 
minor axis. Unlike this behavior, the amplitude of the sec- 
ond harmonic produced by a density wave should weakly 
correlate with the orientation of the disk. Hence, to carry 
out a Fourier analysis of the azimuthal brightness distri- 
bution under initial assumption of zero inclination, and 
to draw the line of maxima of the second harmonic, one 
can get the orientation of the disk major axis (see also 
Grosbol, 1985, who used a similar approach for the outer 
parts of real galaxies). The most clearly delineated orien- 
tation of the line of maxima along the major axis would 
be observed in the regions of a weak density wave (that is 
for a poorly traced spiral structure or a bar), and even if 
the width of each such region is rather small, coincident 
directions of the lines of the harmonic maxima over all 
these regions enable us to estimate an orientation of the 
major axis PAq with good accuracy. It is essential that 
the method does not demand these free-of-distortion re- 
gions to be preliminarily found in a galaxy: they reveal 
themselves in the process of the second harmonic phase 
calculations. 

The inclination angle i may be found (in the absence 
of the non-axisymmetrical structure) as the angle where 
the phase of maxima of the second harmonic of the depro- 
jected image switches from the major axis to the minor 
one. In the presence of a density wave the interference of 
two factors (the error of i and the influence of the spiral 
structure) causes the phase of maxima of the second har- 
monic to be shifted into the position between the major 
(or minor) axis and the observed brightness maximum of 
a spiral arm in the given annulus. In the annuli where 
the spiral arms have the lowest contrast (they may be re- 
vealed during the procedure of PAq determination) the 
phase shift should be the least. Let us note that when 
the probed inclination passes through its correct value, 
the line of phase of the second harmonic turns shortly in 
these regions because of the jump (at 90°) of the position 
of the maximum contribution of the error of i in this har- 
monic. It allows us to estimate i as an angle where the 
phase of the second harmonic best fits the observed spi- 
ral pattern (including the regions of low brightness of the 
arms) . 

The photometric method described above has been 
successfully applied to a sample of double-barred galax- 
ies by Moiseev et al. (2004); the reader may find some 
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Kinematics 


Ring 


average 


Photometry 


Lcda 


RC3 


Name 


Vsys 


i. 


PA," 


range 


j,° 


PA," 


i. 


PA," 




i, 


PA" 


i, 


PA" 


IC 1525 


5008 ± 2 


43 ± 5 


230 i 2 


T <C 90 


44 ± 21 


230 ± 4 


43 ± 3 


210 i 5 




48 


200 


44 


200 


NGC 23 




45 ± 4 


174 ± 2 


r < 37 


48 ± 18 


173 ± 9 


45 ± 3 


170 ± 5 




43 


188 


50 


188 


NGC 615 


1875 ± 3 


71 ± 3 


156 ± 2 


8 < r < 85 


66 ± 18 


153 ± 15 


65 ± 3 


159 ±1 r < 


15" 
























156 ± 2, r > 


20" 


73 


205 


66.5 


155 


NGC 972 


1512 ± 2 


62 ± 4 


325 ± 2 
332 ± 4 


r < 40 
40 < r < 45 


53 ± 23 


320 ± 9 


64 ± 1 


332+i 




65 


332 


59 


332 


NGC 1084 


1438 ± 3 


56 ± 3 


216 ± 2 


r < 60 


54 ± 15 


215 ± 4 


56 ± 2 


217 ± 2 




46 


295 


56 


215 


NGC 1134 


3650 ± 3 


53 ± 4 


328 ± 2 


r < 50 


54 ± 18 


326 ± 5 


42 ±3, r < 40" 


















67 ± 3 


330 ± 3 


50 < r < 90 






67 ± 2, r > 50" 


326 ± 3 




76 


328 


70 


328 


NGC 2964 


1320 ± 4 


47 ± 5 


278 ± 2 


r < 38 


47 ± 24 


276 ± 5 


42 ± 2 


276 ± 2 




58 


277 


57 


277 


NGC 3583 


2068 ± 3 


46 ± 6 


121 ± 3 


r < 70 


45 ± 24 


122 ± 12 


49 ± 2 


124 ± 2 




55 


125 


50 


125 


NGC 3893 


953 ± 3 


42 ± 4 


345 ± 2 


r < 70 


43 ± 18 


345 ± 3 


40 ± 4 


346 ± 3 




32 


345 


52 


345 


NGC 4100 


1150 ± 3 


60 ± 5 
72 ± 2 


358 ± 3 
350 ± 2 


r < 12 
25 < r < 90 


63 ± 17 


356 ± 10 


73 ± 1.5 


353 ± 2, r < 
346 ± 2 


10" 


78 


347 


71 


347 


NGC 4136 


586 ± 2 


25 ± 9 


252 ± 4 


r < 95 


33 ± 16 


252 ± 6 


!8-io 


252 ± 8 




22 


179 


0-29 




NGC 4414 


742 ± 2 


56 ± 4 


155 ± 2 


r < 90 


50 ± 14 


156 ± 2 


56±^ 


157 ± 3 




54 


155 


56 


155 


NGC 4814 


2520 ± 4 


50 ± 5 


297 ± 3 


r < 90 


45 ± 15 


297 ± 6 


51 ± 3 


300 ± 3 




46 


315 


42 


315 


NGC 5371 


2618 ± 5 


47 ± 8 


195 ± 4 


r < 110 


49 ± 21 


193 ± 10 


45 ± 5 


195 ± 3 




48 


188 






NGC 6643 


1487 ± 3 


61 ± 4 


40 ± 3 


r < 70 


54 ± 17 


39 ± 6 


58 ± 2 


37 ± 2 




57.5 


38 


60 


38 



details of practical use of this method in the cited paper. 
The results of "photometric" estimation of the orienta- 
tion parameters for the present sample galaxies are given 
in the fifth column of Table The confidence intervals 
characterize an uncertainty caused by the non-vanishing 
amplitude of the density wave. For the imaging data we 
used (red) continuum brightness distributions obtained 
during our observations with IFP. For three galaxies the 
published I-band images obtained at the lm "Zeiss-1000" 
telescope of the SAO RAS were used. These are NGC 615 
(Sil'chenko et al.HEHJ), NGC 1084 (MoiseevEHni), and 
NGC 1134 (Bizyaev et al. TMity . The R-band image of 
NGC 4136 obtained at the 1.1m telescope of the Lowell 
Observatory was taken from the Digital Atlas of Nearby 
Galaxies (Frei et al. I1996|) . For some galaxies raw photo- 
metric data from the ING Archive of the UK Astronomy 
Data Centre were also retrieved. These are the observa- 
tions of NGC 23 in R-band and NGC 6643 in I-band at 
the lm Jacobus Kapteyn Telescope and the observations 
of NGC 4100 and NGC 6643 in R-band with the 2.5m 
Isaac Newton Telescope. The K-image of NGC 972 was 
obtained with the UKIRT telescope (Hawaii) and is kindly 
provided by Dr. S.D.Ryder. 

3.3. A comparison of "kinematic" and "photometric" 
estimates of orientation parameters 

It is convenient to compare the estimates of i and PAq 
presented in Table |31 as follows: 1) by comparing our 
"kinematic" parameters obtained for the whole velocity 
field with those taken as the average of the parameters 
found for individual rings in a tilted-rings model; 2) by 
comparing our "photometric" parameters with the incli- 
nation and position angles contained in the LEDA and 
RC3 Catalogue, and 3) by comparing the "kinematic" es- 
timates with the "photometric" ones. 



3.3.1. The whole field versus ring average method. 

The confident intervals for the parameters i and PAq ob- 
tained as the means over individual rings are much larger 
than those from our kinematic estimations. Although our 
whole-field estimates formally fall into the wide confidence 
intervals of the "tilted-rings model" results, so that both 
could be considered as "consistent" 4 , the accuracy of the 
parameters obtained for the whole velocity field is signif- 
icantly higher than that using the average of parameters 
found for individual rings in a tilted-rings model 5 . 

3.3.2. Photometric method taking into account spiral 
structure versus outer-isophote methods. 

An agreement between our estimates of i and PAq ob- 
tained by the photometric method taking into account 
the regular spiral structure with those listed in the LEDA 
database and RC3 catalogue appears to be rather mod- 
erate. Note however that the accuracy of the LEDA and 
RC3 data is certainly lower than ours as they are based 
on small-scale Palomar images of galaxies. Concerning 
LEDA, one of the parameters i and PAq are "consistent" 
(within v^Ai) for about a half of the objects, and only for 
3 galaxies (NGC 972, 4414, and 6643) the "consistency" 
exists for both orientation parameters. A comparison with 
the RC3 catalogue gives better results: inclination esti- 
mates are "consistent" for 8 of 14 and PAq - for 10 of 13 



4 Formally speaking, two values a ± Ai and 6± A2 are 'con- 
sistent' if \a-b\ < VAi2 + A 2 2. 

5 We note above that it is not correct to use the statistical 
error of the means as the measure of the confidence interval 
when averaging over individual rings. However, if one does it 
formally, the difference between the angles estimated for the 
whole field and by averaging over rings would be greater than 
the confidence interval. 
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common galaxies. Both parameters are 'consistent' simul- 
taneously for 6 galaxies. 

Nine galaxies of our sample are also included in the 
list of objects for which Grosbol (1985) estimated the 
orientation parameters by analyzing the second Fourier 
harmonic of the azimuthal brightness variations on small- 
scale Palomar images at several radial bins (IC 1525, NGC 
23, 2964, 3583, 3893, 4136, 4414, 4814, 5371). The root 
mean square difference between his values and ours ob- 
tained by the photometric method is 6° for i and 11° for 
PAq. The root mean square differences between Grosbol's 
and LEDA data are slightly higher, 8° for i and 12° for 
PAq (even if to exclude NGC 4136 with strongly dis- 
crepant estimates of PA ). 

3.3.3. Kinematic versus photometric parameters. 

A comparison of the estimates of PAq we obtained by 
the kinematic and photometric methods described in the 
previous subsections reveals discordant results in 3 of 15 
galaxies (IC 1525, NGC 972, and NGC 4100). For the first 
galaxy the disagreement is most likely due to the presence 
of a rather strong bar. For NGC 972 the discordant esti- 
mates may also be related to the bar. Although this galaxy 
is commonly classified as a non-barred one, the presence 
of the inner bar is found in the infrared (see the next sec- 
tion). In NGC 4100 the difference of PA found by two 
methods is rather low (if to exclude the inner part, the 
kinematic and photometric PAqS are 350 ± 2 and 346 ± 2 
degrees respectively) and may be coincidental. 

The non-circular motion caused by spiral density waves 
can distort our kinematic estimates of PAq, but these dis- 
tortions are different in different parts of the disk and 
as a result they partly compensate each other when the 
whole velocity field is taken for the PAq determination. 
The presence of a bar may produce stronger non-circular 
motion with the preferred orientation that are hard to 
compensate by averaging. It may explain the differences 
between the orientation parameter estimates obtained in 
some barred galaxies by the two methods. Figs. and 21 
demonstrate the shapes of the lines of maxima of the sec- 
ond Fourier harmonics of the surface brightness for the 
barred galaxies IC 1525 and NGC 972 calculated using 
different estimates of PAq. As one may see, the shape of 
the lines calculated with the PAo kinematically estimated 
(left) disagrees with the observed spiral arms. Instead, 
these lines tend to oscillate around some preferred direc- 
tion which is close to the orientation of the major axis ac- 
cording to the photometric estimate of PAq. The shapes 
of the lines calculated with the latter PAq values (right) 
are in accordance with the spiral arms traced by eye. All 
these considerations allow us to suggest that the true ori- 
entation parameters of these two barred galaxies should 
be closer to those derived from the photometric data. 

Concerning the inclination i, its values found by us 
from kinematic and photometric data are "consistent" 
in 13 galaxies of 15. Those two galaxies (NGC 615 and 



NGC 1134) where the agreement is absent demonstrate 
rather complex asymmetric structure. 

Indeed, we would advise use of both methods together 
when trying to estimate galaxy orientation parameters: a 
convincing reason is that the non-axisymmetry of a galac- 
tic disk - a bar or a strong density wave - affects the pa- 
rameter estimations by photometric and kinematic meth- 
ods in a different way. 

3.4. Rotation curves in a model of pure circular gas 
motion 

FigEl shows the results of data processing for the galaxy 
IC 1525. Similar illustrations for other galaxies of our sam- 
ple are presented in the electronic version of the article. 
For each galaxy the first three plots show the images in 
the red continuum and Ha (gray-scaled in proportion to 
the square root of the intensity) , and line-of-sight velocity 
field of ionized gas with the isovelocity lines superimposed 
with a step of 50kms _1 (the thick line corresponds to the 
systemic velocity V sys ). The coordinates of the center of 
a galaxy and its systemic velocity were found from kine- 
matic data as the point of the best mirror symmetry of 
the velocity fields for all pairs of pixels taken on opposite 
sides with respect to the center. A small variation of the 
center coordinates has little effect on the orientation pa- 
rameter estimates because they are influenced by different 
Fourier harmonics of the velocity field (see Lyakhovich et 
al. 1997). Additionally, we have determined the positions 
of the center from both photometric and kinematic data 
independently for several galaxies to be sure that they 
coincide within 1-2 pixels and it is always the case. 

In the last plot the rotation curve is shown (the upper 
curve). It was calculated from the observed velocity field 
(see Eq. |2Jl as the best fit curve of circular rotation for 
the "kinematic" orientation angles i and PAq, which were 
assumed to be constant along the radius. The lower curves 
are the variations of i(r) and PAq{t) obtained for the 
individual rings in a tilted-rings model. 

Since each elliptical ring contains from several tens to 
several hundreds of measured points, formal errors of the 
mean V m t(r) in the most cases are lower than 1 — 5kms _1 . 
However, as it was pointed out by Lyakhovich et al. ([1997) 
and Fridman et al. I|1997|l . the rotation curve obtained in 
the frame of the model of pure circular motion may have 
a systematic error of the order of non-circular velocities 
in a density wave. For this reason, with the error bars we 
mark the mean residual velocity in the individual rings 
(\/x 2 ( r i)/Ni, see Eq. 0J. This way is evidently more cor- 
rect than using formally determined errors of V ro t(r) be- 
cause it demonstrates a real confidence in these estimates. 
Genuine circular rotation curves, which may be found by 
methods taking into account the regular non-circular mo- 
tion, will be within the error bars of the curves presented 
here. The large error bars seen at some rotation curves 
give evidence for significant deviations from the circular 
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Fig. 3. A position of maxima of the second Fourier harmonic of the surface brightness calculated for different choice 
of PAq superimposed onto the image of IC 1525. Left (a, c): PAq is chosen according to the kinematic data analysis. 
Right (b, d): this parameter was taken from the photometry analysis (see the text). Images are in the red continuum 
(above) and in Ha line (below). 

Fig. 4. A position of maxima of the second Fourier harmonic of the surface brightness calculated for a different choice 
of PA superimposed onto the K-band image of NGC 972. Left (a): PA is chosen according to the kinematic data 
analysis. Right (b): this parameter was taken from the photometry analysis. 



motion model rather than for the low accuracy of mea- 
surements. 

4. Results for individual objects 

Figures illustrating the results obtained for all galaxies of 
our sample but the first one, IC 1525, are available in the 
electronic version of the paper. 

4.1. IC 1525 

Below we will use PA and R as polar coordinates in the 
sky plane (where PA is the position angle measured from 
the North counterclockwise and R is the distance from the 
galactic center). 

The continuum image (Fig.[S]a) reveals the presence of 
a prominent bar aligned at PA = 48°. The bar ends are 
encircled by the ring of HII regions which has a radius of 
15"- 17" (Fig. Elb). In the central part of the bar, R < 
5" — 7", due to the superposition of strong absorption lines 
of bulge stars, the Ha emission line sinks in absorption, 
so the velocity field (Fig. 0c) has a hole in the center. 

During the observation of this galaxy, the working 
spectral range of IFP contained not only Ha, but also 
the emission line [NII]A6583 from the neighboring inter- 
ference order, the most prominent line in the circumnu- 
clear region (see Fig. Qb). We have tried to estimate the 
line-of-sight velocity of the ionized gas in the central re- 
gion lacking the Ha emission by using the [Nil] line, but 
the velocity field obtained in such a way deviates strongly 
from the simple extrapolation inward of the mean rotation 
field obtained from the Ha measurements further from the 
center. We think that the anomalous velocities obtained 
from the [Nil] measurements in the region r < 5" may re- 
sult from the impossibility to properly take into account 
the effect of the stellar absorption line Ha disturbing the 
profile of [Nil] emission line because of the overlapping of 
interference orders. To treat the gas motion in the center 
of IC 1525 correctly, new observations of emission lines 
free of stellar absorptions are needed. 

The residual velocity distribution, obtained by the sub- 
traction of the mean rotation field from the observed ve- 
locity field, demonstrates several regions where a projec- 
tion of non-circular velocities onto the line of sight ex- 
ceeds 50kms -1 . The sizes of these regions corrected for 
the beam smearing are 1" — 5" (up to 1.6 kpc for the 
adopted distance of 67 Mpc). In the cylindrical coordinate 



frame on the sky plane, at the locus R = 28", PA — 190° 
the residual line-of-sight velocity is 50— SOkms" 1 ; at the 
loci (R = 33", PA = 34°), (R = 22", PA = 105°), 
and (R — 35", PA = 75°) the residuals have the op- 
posite sign: —(40 — 70)kms~ 1 . It is interesting that the 
intensity ratio [Nil] /Ha is enhanced in these regions by 
a factor of 1.5-2 with respect to the nearest outskirts of 
these local velocity anomalies (see Fig. ^a). At the lo- 
cus (R = 26", PA = 254°) the relative intensity of [Nil] is 
quite normal, but the width of the Ha emission line profile 
is FWHM = 200- 220 kms" 1 , (after deconvolution with 
the instrumental contour), and the residual line-of-sight 
velocities are —(50 — 80)kms _1 . All five sites of anoma- 
lous velocities are located beyond bright HII regions and 
are not related spatially to the spiral arms. 

A sharp steepening of the very outer part of the veloc- 
ity curve in this and some other galaxies discussed below 
is evidently an artifact; the result of a small number of re- 
liably measured points (often located near the minor axis 
of a galaxy) and/or of a bad agreement with a model of 
co-planar circular motion. 

4.2. NGC 23 (Mrk 545) 

The rotation curve of the galaxy (Fig. HJd) possesses a lo- 
cal maximum first noted by Afanasiev et al. (1988a). We 
have prolonged this rotation curve to the radius twice as 
large as that reached by Afanasiev et al. (1988a). Since at 
r > 30" there is only one bright "spiral arm" that con- 
tributes to the rotation velocity measurements, the local 
minimum of the rotation curve at r « 40" is a possible ar- 
tifact related to non-circular motion within this arm. The 
inclination of individual rings in a tilted-rings model vary 
strongly along the radius that evidently reflects significant 
deviations from circular motion of about 20 — 30kms , 
nevertheless the inclination i and the position angle of 
line-of-nodes PAq we found for this galaxy (see Table |3J 
remain to be consistent within their error with the ear- 
lier photometric and kinematic estimates (Afanasiev et al. 
1988a). The PAq of the kinematic major axis orientation 
in the inner part of NGC 23 also coincides within a few 
degrees with the estimate of Afanasiev et al. (1991b) made 
with the integral field (multipupil) spectrograph MPFS. 

The problem of possible bar presence in NGC 23 re- 
mains unsolved. According to RC3, the galaxy is classified 
as SB, and de Jong (1996), by applying a 2D decompo- 
sition method to its brightness map, finds a bar with a 
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length of 29" - 32" aligned along PA = 154°. But the 
velocity field obtained by us does not demonstrate any 
noticeable turn of the kinematic major axis (Figs. EJc and 
Eld), and there is also no sign of any isovelocity twist- 
ing when passing from the outer to the inner region of 
the galaxy (Afanasiev et al. 1991b). The local small-scale 
twisting of the isovelocities at r — 5" — 8" (Fig. EJd) is 
rather related to a small (mini- ) bar with a length of 
about 10", which is implied by the CCD- observations of 
Chapelon et al. (1999). 

In the circumnuclear region (r < 5") the emission line 
Ha is broadened up to FWHM = 200-300 kms -1 , which 
was already noted by De Robertis & Shaw (1988) from 
their high-resolution spectroscopy. Their Fig. 1 demon- 
strates two local maxima of FWHM located symmetri- 
cally with respect to the nucleus at r — 6" — 8". Our 
IFP spectra of these spots reveal strong either red or 
blue wings of the Ha line profiles. The Gauss analysis 
of these profiles allows us to derive "non-circular" veloc- 
ity components shifted with respect to the main compo- 
nent by 130 - 200 kms" 1 at (R = 6", PA = 323°) and 
by -(170 - 200) kms" 1 at (R = 3", PA = 190°). The 
profiles of the emission lines in these regions are shown 
in Figs. 0c and d. The estimated diameters of these 
regions are 3" — 4", being close to our spatial resolu- 
tion element. We should note that it is at the location 
of the "blue" component where Afanasiev et al. (1991b) 
had found non-circular motion of the ionized gas with ve- 
locities of 260 — 300 kms -1 , but the asymmetry of the 
emission line profiles was missed by them because of their 
low spectral resolution. The same authors have argued 
that NGC 23 may be at the post-Seyfert stage. In this 
case the anomalous gas motion may be related to the ac- 
tivity, or to a starburst in the nucleus. But the emission 
line intensity ratios in the nucleus of NGC 23 evidence 
rather in favor of HH-region character of excitation and 
are not consistent with the hypothesis of Sy2 or LINER 
(De Robertis & Shaw 1988; Contini et al. 1998). Since the 
sites of the anomalous velocities are located at the ends of 
a nuclear minibar (see above), the presence of the second 
velocity component may be related to non-circular motion 
in the bar. But in this case one needs to explain a very 
high degree of non-circular motion - more than 100 kms -1 
in projection along the line of sight. 

4.3. NGC 615 

In the circumnuclear region (r < 10" — 15") the broad 
absorption line Ha is superposed on the emission line 
[NII]A6548 from the neighboring interference order that 
prevents an exact determination of the line-of-sight veloc- 
ities even by a Gaussian analysis. We start our measure- 
ments beyond the radius of the circumnuclear local rota- 
tion velocity maximum (at r « 5") detected by Afanasiev 
et al. (1988a). Beyond r = 10" the relative intensity of 
the nitrogen emission lines falls, so they do not affect the 
profile of Ha any more. 



The most prominent feature of the monochromatic Ha 
image (Fig. 0b) is a pair of HII regions located sym- 
metrically with respect to the nucleus at the radius of 
20" — 23". They are slightly elongated along the radius 
and look like shock fronts at the edges of a bar (though 
the galaxy is classified as unbarred) . At the same distance 
from the center the gas velocity field is disturbed by non- 
circular motion resulting in a turn of the kinematic major 
axis by more than 20° at r = 15" — 30". We have found 
here several spots with the line-of-sight velocities differ- 
ing by 40 — 50 kms -1 from the mean rotation field. These 
kinematic details may be evidence for a triaxial potential 
within r = 25" — 30" in this galaxy. Indeed, the recent pho- 
tometric analysis of Sil'chenko et al. (2001) has revealed 
that NGC 615 possesses a separate inner compact disk of 
oval shape especially notable at r = 15" — 30" where it 
introduces non-monotonic changes of isophote orientation 
and ellipticity. 

4.4. NGC 972 

The near-Ha continuum image (Fig. |H|a) lacks a spiral 
structure, but on the Ha image (Fig. 0b) the HII re- 
gions concentrate toward a pseudo-ring with the radius 
of about 30". Several works which analyze CCD-data on 
NGC 972 (Ravindranath & Prabnu 1998, Zasov & Moiseev 
1998) also contain remarks about the absence of a no- 
ticeable spiral structure on the broad-band optical im- 
ages and about a large amount of dust within the ring 
of HII regions. However, near- infrared images obtained 
through JHK-Glters reveal a prominent two-armed spi- 
ral pattern and a small bar with the size of 10" elongated 
in PA = 175° - 180° (Mayya et al. 1998; Zasov & Moiseev 
1998). This peculiarity makes the morphological classifica- 
tion of the galaxy very ambiguous; e.g. Mayya et al. (1998) 
classify it as SABdpec while according to NED (RC3) it 
is Sab. 

The isovelocities of the ionized-gas velocity field within 
the central 10" have a characteristic S-shape (Fig. [He), 
evidently induced by a presence of the bar seen in NIR. 
A turn of the kinematic major axis in the center exceeds 
30° (Fig. Eld); the non-circular (radial?) gas streams in the 
bar are so significant that a cross-section of the velocity 
field in PA « 30° gives an impression of counterrotation. 
We should note here that the presence of a minibar in 
the center of NGC 972 had been first claimed by Zasov &; 
Sil'chenko (1996) from kinematic arguments. 

The orientation parameters found by us for the outer 
part of the disk (Fig. Eld) are consistent within the errors 
with both photometric (Ravindranath & Prabnu 1998, 
Zasov & Moiseev 1998) and kinematic results from our 
earlier long-slit observations (see Zasov & Moiseev 1998, 
who used the data from Afanasiev et al. 1991a). The ro- 
tation curve derived in this work also agrees rather well 
with that obtained by Zasov & Moiseev (1998). 

Ravindranath & Prabnu (1998) noted a high rate of 
star formation in the central part of NGC 972, comparable 
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even to that of the well-known starburst galaxies M 82 and 
NGC 253. We see indirect consequences of this intense star 
formation when analyzing our Fabry-Perot spectra: many 
HII regions stand out by their negative residual veloci- 
ties of about —(15 — 20) kms -1 , which is probably related 
to expanding bubbles of ionized gas. A negative sign of 
residuals may be accounted for by the presence of inter- 
nal absorption, dimming the radiation of receding parts 
of the bubbles. In the nucleus, at r < 2" — 3", which is 
close to our spatial resolution, the Ha emission-line profile 
has mostly a two-component structure; it can be related 
either to gas outflow from the starburst nucleus or to non- 
circular gas motion at the edges of the minibar (analogous 
to NGC 1084, see below). 

4.5. NGC 1084 

The emission line [NII]A6583 from the neighboring inter- 
ference order is seen throughout the galaxy; it has allowed 
us to study the kinematics of the ionized gas by using two 
emission lines separately instead of a single Ha. From the 
earlier observations by Afanasiev et al. (1988a), two re- 
gions with strong velocity deviations from the mean rota- 
tion were found. Particularly in the central region, r < 5", 
the rotation velocities obtained by measuring the forbid- 
den emission lines [Nil] and [SII] were systematically lower 
than those obtained by measuring Ha. The authors had 
interpreted this peculiarity as evidence for minibar pres- 
ence in the center of NGC 1084. Besides that, non-circular 
motions were found at the periphery of the galaxy, at 
r = 40" — 50" to the north-east from the center, near 
an elongated arc-like spur outgoing from the spiral arm 
near a giant superassociation (Fig. 0b) . 

Our IFP data confirms the existence of these kinemat- 
ically distinct regions. Near the spur, the profile of the 
Ha emission line is mostly two-component: one compo- 
nent traces the mean rotation and the other is shifted by 
±100 — ±150 kms -1 . The appearance of the non-circular 
component is supplemented by a sharp increase of the 
[Nil] /Ha ratio by a factor of 2 — 6 with respect to the 
normal one in the neighboring HII regions. The areas of 
non-circular motion have typical sizes of 1-2 kpc and are 
located between the HII regions. A possible interpreta- 
tions of these non-circular motions of ionized gas near the 
spur is discussed elsewhere (Moiseev 2000). 

In the circumnuclear region, r < 5", the profiles of 
the Ha and [Nil] emission lines also demonstrate a two- 
component structure, with the velocity separation of the 
components about 50 kms -1 ( Figs. ^ e and f). The 
Gaussian analysis of the Ha profile has allowed us to trace 
azimuthal variations of the observed velocity gradients for 
each component. We have found that the brighter compo- 
nent traces the circular rotation of the ionized gas with 
the kinematic major axis close to the major axis of the 
galaxy, whereas the fainter component of both lines cor- 
responds to the gas radial motion because its projected 
velocity gradient is zero near the line of nodes. 



In the distance range of 5"-20" on the opposite 
sides of the nucleus, in the direction PA « 60°, the ex- 
tended linear front is revealed. At both sides of this front 
the character of asymmetry of the Ha and [Nil] profiles 
changes abruptly: the red wings are replaced by blue ones. 
The line-of-sight residual velocities of the main line com- 
ponents change their sign across the front jumping by 
80 - 100 kms" 1 . The width of the front is 3" - 4", so 
it remains practically unresolved. Such a combination of 
features evidences that we are dealing with supersonic gas 
flow in a strong shock front. 

Features similar to those observed in NGC 1084 (radial 
gas flows, narrow shock fronts) are usually treated as a 
result of bar influence; but neither our continuum image of 
NGC 1084 (Fig. 0a), nor broad-band photometric images 
of the galaxy (Zasov & Moiseev 1999, Moiseev 2000) reveal 
any signatures of a bar. We should also note that the turn 
of the kinematic major axis in the radius range of the 
shock fronts (5" - 20") is rather small (Fig. 0d), which 
implies the absence of a high contrast triaxial potential. 

A whole complex of peculiarities observed in the center 
of NGC 1084 may be interpreted as a result of a weak 
triaxial potential (a small triaxial bulge). Similar shock 
fronts at the edges of a slightly triaxial bulge have been 
also observed in Sb galaxy NGC 2841 by Afanasiev & 
Sil'chenko (1999), but in that galaxy the bulge is much 
more massive. 

4.6. NGC 1134 

Due to its complex inner structure (multiple dust fil- 
aments, etc.) the galaxy is included in Arp's Atlas of 
Peculiar Galaxies (1966), though it is not an interacting 
system. In the area r < 20" — 30", HII regions trace a 
three-arm spiral (Fig. 1101b). but in the outer part of the 
galaxy the spiral pattern is two-armed. The three-arm spi- 
ral structure in the center of NGC 1134 can also be seen 
when analyzing its direct images by constructing maps of 
a color Q-parameter free of dust influence (Bizyaev et al. 
2001). 

The mean orientation of the line of nodes of the 
disk (Fig. HOkO is consistent with the earlier photometric 
and kinematic estimates (Afanasiev et al. 1991b), though 
our estimate of the inclination is significantly lower than 
i = 72° used by these authors. At small radii, r < 12", the 
kinematic major axis (Fig. II Ol d) turns abruptly by more 
than 15°, so that the isovelocities fFig. 1101c) demonstrate 
an S-like shape typical for a bar potential. Above all, on 
both sides of the nucleus, at r = 4" — 7" there are regions 
where Ha emission-line profiles are two-component, with 
a separation between the components of 150 — 200 kms" 1 
( see Figs. g and ED h. These regions of the two- 
component Ha emission line are elongated along the di- 
rection of PA as 20° — 30° and resemble the shock edges 
in NGC 1084. Close to these shock fronts, at the minor 
axis of the galaxy, one can note two spots where the sec- 
ondary velocity component disguises itself as an asymmet- 
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ric wing of the Ha emission-line profiles, and the velocities 
obtained for the peak of profiles differ from the mean rota- 
tion velocities by ±(80—100) kms" 1 . Since the anomalous 
velocities are observed near the minor axis of the galaxy, 
they may be explained as the presence of radial gas flows. 
We should mention that for the first time the regions of 
the non-circular gas motion on both sides of the nucleus 
in NGC 1134 were detected by Afanasiev et al. (1991b). 

All the kinematic peculiarities of the central region 
of NGC 1134 mentioned above can be explained by the 
hypothesis of a bar aligned with the minor axis of the 
galaxy. 

The mean rotation curve calculated by us (Fig. HTHd'1 
agrees with the one reported earlier by Afanasiev et al. 
(1991b): the solid-body area up to r = 7", a low-contrast 
peak, and after that a velocity plateau. At r = 65" — 75" 
and at r = 90" — 110" there are two emission "islands" 
in the disk, emerging by 70 — 90kms~ 1 from the extrap- 
olated flat velocity of rotation, which are responsible for 
the sharp jumps in the calculated velocity curve at large 
distances from the center. These anomalous details can 
be related either to a gas-disk orientation change at large 
radii (a "warp"), or to the possibility that these HII re- 
gions do not belong to the global disk plane, being dwarf 
satellites of NGC 1134. To clarify these possibilities, the 
results of HI observations of the outer disk kinematics are 
needed. 

4. 7. NGC 2964 

The galaxy possesses a well-developed two-arm spiral 
structure; HII regions at r = 30" — 35" are concentrated 
to a pseudo-ring fFigs. ITTla and II lib). It has a close com- 
panion, NGC 2968, at the distance of 5.'8, but there are 
no evident signatures of interaction of galaxies. 

At r = 1" — 5" the kinematic major axis (Fig. II lid) 
turns abruptly by more than 15°, which can be assumed 
to be the effect of a small bar (note that there is a remark 
about the "very faint bar" in RC1). The high-resolution 
image of the central 18" x 18" region of NGC 2964 obtained 
by HST and published by Carollo et al. (1997) reveals the 
presence of a complex pattern of dust filaments tracing 
a two-arm spiral and penetrating into the center of the 
galaxy. Similar dust spirals were found by HST imagers 
in the centers of many galaxies during the last years (e.g. 
Martini & Pogge 1999, Regan & Mulchaey 1999, Barth et 
al. 1995). 

There is a region between the nucleus and a spiral arm, 
at r = 5" — 12" to the south of the nucleus, where Ha 
emission-line profiles are clearly two-component ( Fig.^i)- 
The main line component traces the mean rotation; the 
secondary component is redshifted by ~ 150 kms -1 , some- 
times being as strong as the main one. At the outskirts 
of this region the Ha profiles often possess a red wing, so 
at the residual velocity map the whole area located be- 
tween bright HII regions is distinguished by strong non- 
circular motion reaching 50 — 70 kms . The diameter of 



the anomalous area is 10"— 15" (1 — 1.5 kpc for the adopted 
galaxy distance of 18 Mpc). As this area is located at the 
galactic minor axis, the velocity residuals will be caused 
by radial gas motions if they are confined to the plane of 
the galactic disk. 

4.8. NGC 3583 

This galaxy is a member of the pair of galaxies shared 
with NGC 3577 (their separation is 5'), but there is no 
clear evidence of interaction. A small elliptical companion 
may be found at a distance of about 1' (see the note in 
RC2 catalog). The continuum image of the galaxy reveals 
the presence of a bar which can be traced up to r w 30" 
and is aligned with PA = 80°; it also contains a lot of HII 
regions fFigs. IT2la and 1121b). Isovelocities close to the bar 
have a typical S-like shape (Fig. 1121c) , the kinematic ma- 
jor axis turns here by more than 20°, which is evidence 
of a significant contribution of non-circular motion into 
the line-of-sight velocity field (Fig. I12ld). Indeed, within 
the bar, in the radius range of 2" — 18", the line-of-sight 
velocities differ from the mean rotation velocity field by 
±(50 — 80) kms -1 . The 'blue" and "red" residual veloc- 
ity spots are located symmetrically on both sides of the 
center, almost along the minor axis of the galactic disk, 
which enables us to propose a radial gas flow due to the 
bar. The presence of a broad stellar Ha absorption line in 
the circumnuclear region of this galaxy introduces some 
uncertainty into the line-of-sight velocity estimates, but in 
comparison to other analogous cases (IC 1525, NGC 4414, 
etc), this effect is not as severe. 

4.9. NGC 3893 

The galaxy forms a pair with NGC 3896 without any sign 
of tidal interaction (their separation is 3.'9). The rotation 
curve rises all along the radius, having no clear maxi- 
mum (Fig. I13l d). The turn of the kinematic major axis at 
r < 20" is related to the velocity field perturbation caused 
by a small bar aligned at PA » 20°, which is barely seen 
in the continuum image (Fig. 1131a) . The variations of PAq 
at r = 20" — 50" may be a spiral arm effect. At the radii 
of 37" - 45" in the direction of PA = 140°, at the inner 
edge of the spiral arm, the Ha emission-line profiles have 
a double-peaked structure ( see Fig. Q j), the secondary 
(non-circular) component being redshifted by ±50 kms -1 . 
The region where the secondary velocity component ap- 
pears has a diameter of 8" — 9" (about 500 pc for the 
adopted distance of 13 Mpc), is located between bright 
HII regions, and is related to local dynamical processes in 
the gaseous disk rather than to the spiral arms. 

The orientation angles of this galaxy were also esti- 
mated from the Ha velocity field by Garrido et al. (2002), 
by minimizing the dispersion of the points along the rota- 
tion curve. The value found by them, PAq = 347°, is close 
to the 345° ± 2° accepted in this paper, but the difference 
of inclination angles (30° versus 45° ± 5°) is significant. 
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Note however that the accuracy of measurements of i given 
by Garrido et al. (2002) is rather low, ±(5° - 10°). 

4.10. NGC4100 

This galaxy is a member of the Ursa Major cluster 
(Odenwald 1986), although it has no noticeable neighbors 
- the distance to the nearest luminous cluster galaxy is no 
less than 210 kpc (see the map of the cluster in Tully et al. 

1996) . Both CCD images and surface brightness profiles 
reveal a central "lens" between the radii of 15" and 70" 
containing a pair of tightly wound spiral arms. Beyond 
the "lens" , the surface brightness falls abruptly, although 
some interesting details - particularly two low-contrast 
symmetric arms - can still be seen in the outer parts of 
the disk. 

Our IFP data covers only the region of the central 
"lens"; here the Ha emission is confined to the spiral 
arms and to the very bright circumnuclear region, r < 12" 
fFig.lTflb). According to Pogge (1987), this circumnuclear 
region concentrates more than 20% of the total Ha emis- 
sion of the galaxy. The emission-line intensity ratios in 
the nuclear spectrum are typical for HII regions (Ho et al. 

1997) , hence a high Ha brightness of the nucleus may be 
attributed to intense star formation burst, distinguished 
in the background of the more quiescent galactic disk. 

A kinematic study of this galaxy's ionized gas with a 
long-slit spectrograph (Afanasiev et al. 1988b, 1992) has 
shown that the circumnuclear region of NGC 4100 is dy- 
namically decoupled from the rest of the disk: its rotation 
curve demonstrates two local extrema at r = ±5", with 
a rotation velocity gradient of about 300 kms _1 kpc _1 be- 
tween them, which may be evidence of a highly concen- 
trated bulge (Afanasiev et al. 1988b). Above all, Afanasiev 
et al. (1992) noted that the kinematic major axis in the 
central region, r < 5" — 6", was turned by about 25° 
with respect to the outer isophote orientation. The au- 
thors treated this fact as a result of minibar influence in 
the center of this galaxy. 

Our analysis of the velocity field of t he ionized gas in 
NGC 4100 (Fig. 1141c) confirms the dynamical decoupling 
of the circumnuclear region, although the local maximum 
of the rotation curve (Fig. I14H1 is not as prominent as it 
appeared followed from Afanasiev et al. data. Besides the 
turn of the kinematic major axis, from PAq — 350°±2° for 
the outer galaxy to 358° ±3° for the inner one fFig.H4ld). 
the inclination determined from the kinematics in the 
tilted-rings model changes abruptly from i = 72° ± 2° 
for the outer galaxy to i = 60° ± 5° for the inner one. As 
it was noted in the third section of this paper, the change 
of i in a tilted-rings model may be just a reflection of the 
ordered non-circular coplanar motion of a gas. In general, 
presence of a bar may easily imitate the inclination varia- 
tions in a line-of-sight velocity field (Moiseev & Mustsevoi 
2000). But in this particular galaxy, the photometric ma- 
jor axis turns in the same direction as the kinematic one, 
so we cannot exclude the presence of an inclined gas disk in 



the circumnuclear region of NGC 4100 (Zasov & Moiseev 
1999). The radius of this disk is less than 11" - 12" (0.8 
kpc for the adopted distance of 16 Mpc). Depending on 
which side of the circumnuclear disk is nearest to us, the 
inclination of the central disk to the global galactic plane 
must be either 25° or 55°. 

Our rotation curve for NGC 4100 (Fig. I14ld) demon- 
strates a monotonic rise of the rotation velocity up to 
the spatial limits of the measurements, which is consistent 
with HI observations of Verheijen (1996), demonstrating 
that the rotation velocity rises up to r = 90"; in the very 
outer part of the galaxy the rotation velocity reaches a 
plateau. 

4.11. NGC 4136 

According to Allsopp (1979), the galaxy is a possible mem- 
ber of the CVnII group. In the central part (r < 15") there 
is a small bar weakly recognizable in our continuum im- 
age (Fig. 1151a). but clearly seen on the CCD images given 
in The Digital Atlas of Nearby Galaxies by Frei et al. 
(1996). The galaxy is seen almost face-on, and as a result, 
the parameters i and PAq are poorly determined (see the 
dispersion of the estimates in Table|3J). For this reason, the 
absolute value of the rotation velocity is rather uncertain. 
The "troughs" at the rotation velocity (Fig. I15ld) may be 
just artifacts. Formally, this uncertainty is expressed by 
large error bars in Fig. ^j] Under such a disk orientation, 
the line-of-sight velocity deviations from the mean rota- 
tion in most cases may be related mostly to gas motion 
along z-coordinate direction. 

At the radii of 4" — 15" the Ha emission is either negli- 
gible, or strongly contaminated by the absorption line due 
to bulge stars; this effect cannot be properly taken into ac- 
count because of the IFP interference order overlapping. 

4.12. NGC 4414 

According to Elmegreen & Elmegreen (1987), NGC 4414 is 
a flocculent galaxy; our continuum image (Fig. 1161 a) does 
not show any traces of a spiral pattern either. However, ob- 
servations in the NIR bands reveal the presence of several 
extended segments of spiral arms populated by old stars 
(Thornley 1996). In the circumnuclear region of NGC 4414 
we see a narrow Ha emission line inside the broad ab- 
sorption produced by a contribution of bulge stars; within 
several arcseconds from the galactic center only the [Nil] 
emission line from the neighboring order of interference is 
seen (Fig. ^ k). But unlike other similar cases (IC 1525, 
NGC 615, etc), for this galaxy we have been able to con- 
struct a velocity field of the ionized gas in the central re- 
gion by subtracting the absorption line superposed and by 
applying a Gaussian analysis to the Ha and [Nil] emission- 
line profiles. As a result, only a small area with a size of 
7" to the south-east of the nucleus remained uncovered by 
our measurements (Fig. 1161 c) . Our orientation parameters 
for NGC 4414 (PA = 156°, i = 56°) are very close to those 



16 



Fridman et al.: The orientation parameters of spiral galaxies. 



obtained by Thornley & Mundy (1997) from the analysis 
of the HI and CO velocity fields (PA = 160°, i = 60°). 
The HI velocity field measured in this work extends be- 
yond three optical radii of the galaxy and demonstrates a 
strong asymmetry of the outer disk at r > 200", although 
no optical satellite is found within the nearest outskirts of 
NGC 4414 (Thornley & Mundy 1997). 

The turn of the kinematic major axis and large error 
bars of the mean rotation velocity estimates at r < 9" 
(Fig. I16k may be caused by non-circular gas motion in 
the minibar area. The presence of a minibar in NGC 4414 
was first suspected by Thornley & Mundy (1997) from 
isovelocity twisting seen in the HI velocity field, though 
the spatial resolution of the HI observations was much 
worse than our resolution of the ionized gas observations. 

4.13. NGC 4814 

The main fraction of the Ha emission in this galaxy is 
confined to the inner disk having a radius of 50"; beyond 
this radius the emission regions are concentrated in a pair 
of spiral arms (Fig. 13b). A mean rotation curve (Fig.ll7ld) 
agrees well with the rotation curve of the inner region, r < 
30", we obtained from the earlier long-slit observations at 
the 6m telescope (Zasov & Sil'chenko 1987; Afanasiev et 
al. 1988b). A turn of the kinematic major axis at r < 15" 
may be related to a small bar, but it is not accompanied 
by the significant continuum isophote turn. Afanasiev et 
al. (1988b) reported non-circular gas motion to the west 
of the nucleus, at r « 20". After subtracting the mean 
rotation field from our velocity data, we have found that 
the residual velocities in this region are about —40 km s _1 , 
evidently produced by an effect of the spiral pattern on 
the kinematics of the gaseous disk. 

4.14. NGC 5371 

This galaxy possesses a modest bar of about 20" — 30" 
length aligned with the minor axis of the galactic disk. 
The bar is barely seen on our continuum image fFig. 1181 a). 
but reveals itself rather well on the other optical images 
(e.g. Gonzalez et al. 1996). Unfortunately the IFP obser- 
vations of this galaxy were made under poor seeing con- 
ditions during a non-photometric night. Due to this, we 
present velocity measurements only for the brightest emis- 
sion regions in the spiral arms fFig. 118( c). In the circum- 
nuclear region the Ha emission line is strongly broadened, 
and the intense [Nil] emission line from the neighboring 
interference order is overlaid (the [Nil] emission is abnor- 
mally strong because NGC 5371 is classified as a LINER 
(Rush et al. 1993)). We have not succeeded in separating 
the contributions of the overlapping broad emission lines 
into the integrated profile. Perhaps due to this uncertainty 
our measurements of the line-of-sight gas velocities in the 
circumnuclear region do not reveal the local rotation ve- 
locity maximum at r = 10" — 15" reported by Zasov & 
Sil'chenko (1987). It is worth noting that the inclination 



angle % we obtained from the velocity field, i = 47° ± 8°, 
is somewhat larger than i — 34°, adopted by Wevers et al. 

(1986) . 

4.15. NGC 6643 

According to the classification of Elmegreen & Elmegreen 

(1987) , the galaxy has a 5th arm class, and therefore is 
almost lacking an ordered spiral structure. Indeed, both 
our continuum and Ha images fFig. 1191 a and !19l b) present 
a flocculent-like galaxy. However, recent NIR observations 
by Elmegreen et al. (1999) reveal the presence of a two- 
arm spiral pattern in NGC 6643 related to an old stellar 
population and implying a dynamical origin. Afanasiev et 
al. (1988b) noted large-scale non-circular gas motion in 
the central part of the galactic disk of this galaxy. Our 
analysis of the velocity field shows that line-of-sight ve- 
locity deviations from the pure rotation field do exist, but 
they do not exceed 40 — SOkms" 1 over the whole disk 
of the galaxy. These motions are probably related to the 
spiral structure; but an exact diagnosis will be possible 
only after a detailed Fourier analysis of the velocity field, 
which is under way. However there are a few regions in the 
galaxy where the Ha emission-line profile cannot be fitted 
by a single Gaussian function. For example, at r = 12" to 
the north-west of the nucleus, between bright HII regions, 
the Ha profiles demonstrate a broad red wing that may be 
evidence of the presence of the secondary "non-circular" 
velocity component of 100 — 200 km s -1 in projection on 
the line of sight (see Fig. ^ X)- The si ze of this anoma- 
lous velocity area is 5" x 10" (about 1 kpc, under the 
adopted distance of 19 Mpc). The effect of this local ve- 
locity anomaly on the whole velocity field is rather strong: 
it is responsible for a turn of the kinematic major axis ob- 
served in the radius range of r = 3" — 12" (Fig. I19k0. If 
we mask this region when analyzing the velocity field, the 
position angle of the kinematic major axis in the central 
region becomes about 40° , in good agreement with the ori- 
entation of the line of nodes of the outer disk (Table |3J) . 
We should also note that the inclination i w 60° found in 
this work is less than i = 68°, used by Afanasiev et al. 
(1988b) but agrees with i = 59° of Bottinelli et al. (1984) 
obtained from photometric data. 

5. Discussion 

Our sample of galaxies is rather widely representative, be- 
cause it includes objects with a variety of kinematic prop- 
erties of gas and with a wide range of structural features. 
The rotation curves in all the cases rise steeply in the in- 
ner parts within r w 10" — 20" (in most cases inside 1 - 
3 kpc from the center) due to the combined potentials of 
the bulge and the disk; further from the center the shape 
of curves looks different: the curves have either a plateau 
with some local velocity variations, or pass local maxima 
(as in the cases of NGC 23, 615, 1134, 3583, 5371). The 
curves of rotation in IC 1525, NGC 972, 1084, 3893, 4100, 
4136, 6643 rise to the limit of the observed part of the 
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galaxies (which, being less than the optical sizes of the 
galaxies, nevertheless contains the brightest part of the 
spiral structure). The opposite tendency is observed in 
NGC 2964; its rotation curve falls beyond the region of 
bright spiral arms. 

The observed large-scale non-circular motions are usu- 
ally connected with a spiral density wave or a bar. In par- 
ticular, these motions result in the artificial PA® and i 
radial variations in a tilted-rings model, even when their 
real values remain constant. Systematic perturbations of 
the line-of-sight velocities in the region of spiral arms are 
most clearly seen in NGC 23, 4814, and 6643 (see, respec- 
tively, Figs.l6ld.H7ld and!19k0. although in other galaxies 
they also appear. 

In the presence of a bar (which can be considered 
another manifestation of a density wave) or some other 
type of triaxial structure (triaxial bulge, oval lense, etc) 
gas clouds follow elongated orbits which results in an ob- 
served turn of the kinematic major axis and leads to a 
well-known S-shaped distortion of isovelocities (see e.g. 
Chevalier & Furenlid 1978). Besides that, in such cases 
a secondary component of the emission line is often ob- 
served; as a rule, it traces radial gas flows along the bar. 
Indeed, these peculiarities are observed in our galaxies, 
but curiously, the velocity anomalies correlate weakly with 
the relative size and contrast of a bar. There are three 
galaxies in our sample classified as SB in both NED and 
LEDA catalogs: IC 1525, NGC 23, and NGC 3583. In the 
former two we observe only a slight hint of bar-related 
anomalies, and only in NGC 3583 the expected distortion 
is clearly observed. In general, among the galaxies with 
the kinematic signatures of a bar there are the galaxies 
with the known photometric bars (SAB and SB, accord- 
ing to NED: IC 1525, NGC 23, NGC 2964, NGC 3583) as 
well as the galaxies previously classified as SA (NGC 615, 
NGC 972, NGC 4414, NGC 4814). However it is well 
known that low-contrast optical, or only infrared, triaxial 
structures may exist in galaxies classified as SA (NGC 972, 
NGC 615). Indeed, in NGC 972 the bar has been detected 
from the NIR surface photometry (Mayya et al. 1998), 
and in NGC 615 an oval inner disk has been found from 
the optical photometric data (Sil'chenko et al. 2001). We 
think also that the local regions with large residual gas 
velocities, situated symmetrically on opposite sides of the 
center and characterized by an enhanced [Nil] emission, 
which are found in NGC 1084 and NGC 1134, may be con- 
nected to the front of a bar-induced shock wave, although 
optical bars are not clearly seen here. 

Besides the large-scale deviations from circular rota- 
tion related to a bar or spiral arms, significant small-scale 
velocity anomalies are also observed in all galaxies we 
studied. They may be divided into three categories. 

The first type refers to the non-circular motions of dif- 
ferent natures in circumnuclear regions (within ten arc- 
seconds from the center). Wide wings or significant shifts 
of spectral lines may be evidence for radial gas motion as 
a result of nuclear activity or an active star formation in 
the nuclei (IC 1525, NGC 23, NGC 972, NGC 1084). One 



should bear in mind, however, that similar manifestations 
(especially from observations with a moderate spatial res- 
olution) may be caused by nuclear mini-bars (NGC 615, 
NGC 972, NGC 2964, NGC 4414) or by an inclined nu- 
clear disk (NGC 4100). More detailed analysis is needed 
to distinguish between these cases (see discussion in Zasov 
& Moiseev 1999). 

The second type of local non-circular motion is ob- 
served in the bright HII regions of a disk. Among our 
galaxies such motions are definitely observed only in the 
most actively star- forming galaxy NGC 972. The velocity 
distortions are not too high (about 20kms~ 1 ). This type 
of motion is most certainly caused by gas outflow as a re- 
sult of heating and blowing by massive stars in the region 
of a local star formation burst. More energetic motions of 
this type are usually observed in star-forming dwarf galax- 
ies (radial expansion with amplitudes of ten to a hundred 
kms -1 , see for example Pustilnik et al. 2001). 

The third type of local non-circular motion is the most 
enigmatic. It takes place far from the nuclei, in the regions 
with a typical diameter of about 0.3-1.5 kpc which exceeds 
the size of bright regions usually covered by a single burst 
of star formation. Abnormally moving gas is revealed by 
its Doppler velocities, which differ from the velocities of 
the surrounding disk or from the velocities expected for 
circular rotation at a given point. The velocity differences 
reach some tens of kms -1 , but sometimes they exceed 
100 kms" 1 . Such regions are found in IC 1525, NGC 1084, 
NGC 2964, NGC 3893, and NGC 6643. These velocity 
anomalies are often accompanied by an appearance of sec- 
ondary emission-line components which indicate the exis- 
tence of at least two independent subsystems of ionized- 
gas clouds possessing normal and abnormal motions on the 
line of sight. Sometimes the presence of a "non-circular" 
velocity component is seen in the wide asymmetric wings 
of the emission line profiles. In IC 1525, NGC 1084, and 
partly in NGC 1134, the non-circular motions are accom- 
panied by an increase in the [Nil] /Ha ratio that argues 
in favor of a shock excitation mechanism (see the detailed 
discussion in Moiseev 2000). 

The regions of these abnormal velocities are usually 
observed off the sites of active star formation, although 
the bright HII regions are often seen near them. It seems 
as if the gas with abnormal velocity is usually located be- 
tween the complexes of HII. Nevertheless it does not mean 
that the "abnormal" gas definitely flows around normally 
rotating islands of star formation: indeed, the extended 
regions of this gas may cover bright HII regions in projec- 
tion, which makes the emission of abnormally moving gas 
non-observable against the bright background. 

These regions may be related to powerful "galactic 
fountains" , where gas thrown out of star formation sites 
along the z coordinate returns, falling toward the galactic 
disk (sec for example Breitschwerdt & Komossa 2000, de 
Alvis 2000). We observe this gas more easily when it is 
projected onto the areas of low gas emissivity. The other 
event which may account for the observed gas motion is 
the interaction of low density gas of a galactic disk with 
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high velocity clouds accreting onto the galaxy. In this case, 
the appearance of star formation sites in this region is a 
consequence, not the cause, of the abnormal velocities of 
the gas. Further studies may distinguish between these 
possibilities. 

6. Conclusions 

The main results of this work may be summarized as fol- 
lows: 

— The presence of regular noncircular motion of gas re- 
lated to spiral density waves may introduce systematic 
errors in the evaluation of the orientation parameters 
i and PA of galactic disks. In particular, it affects the 
commonly used tilted-rings model. The systematic er- 
rors remain even if the parameters found for individual 
rings in this model are averaged. To avoid this effect, 
we force the orientation parameters to be constant all 
over the disk within some radius. Still, the tilted-rings 
model remains a useful diagnostic method of revealing 
non-coplanar or strongly non-axisymmctric gas motion 
at some distance from the center (in circumnuclear re- 
gions or the disk outskirts). 

— For the sample of fifteen observed galaxies the param- 
eters i and PA were estimated using both the kine- 
matic and photometric methods we describe in this 
work; the latter is based on the analysis of the sec- 
ond Fourier harmonic of azimuthal brightness distribu- 
tion in the sky plane. Both kinematic and photometric 
methods of disk orientation parameter determination 
may have their limitations, especially in galaxies with 
strong triaxial (oval or bar) distortions. However, the 
close agreement of the orientation parameters found 
by two independent methods for most of our sam- 
ple galaxies shows that the effect of systematic non- 
circular gas motion does not significantly alter the re- 
sulting parameters obtained from kinematic data for a 
pure circular rotation model. 

— We have obtained the curves of rotation for 15 spiral 
galaxies, assuming the gas motion to be circular, and 
found the confidence intervals for rotation velocities 
at a given radii which characterize possible errors due 
to ignoring the regular non-circular motion related to 
spiral arms and/or a bar. 

— Practically all the galaxies possess local non-circular 
gas motion, which appear either as residual veloci- 
ties left after subtraction of circular velocities from the 
observed velocity fields, or as the complex profiles of 
emission lines, which may be decomposed into circular 
and non-circular components. Local anomalies of the 
velocity fields on a scale of several hundred pc (up to 
1-1.5 kpc) are related to star formation or some other 
mechanisms, which are responsible for the non-circular 
velocities often exceeding 100 km s -1 (see the discus- 
sion in the previous section.) 

— The observed regular large-scale non-circular motion 
in the galaxies have another nature: they are related 



to bars and/or spiral arms of galaxies. In principle, the 
detailed analysis of these density-wave-induced mo- 
tions may allow to restore a full 3D vector velocity 
field of the gas in galactic disks (see the Introduction) . 
This will be done as a future work. 
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7. Appendix. Images and gas velocities of the 
individual galaxies. 

See the text for details 
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Fig. 5. IFP data for galaxy IC 1525. Top - Image in continuum near Ha (left) and monochromatic image in Ha (right), 
a square-root scale is used. Bottom left - velocity field, the thickest black contour corresponds to the systemic velocity, 
other contours are ±50, ± 100, ± 150 ... kms' 1 . A cross marks the dynamical centre. Bottom right - rotation curve 
(Vrot), kinematic position angles (PA), and the inclination (i) (both are obtained in a tilted-rings model). Dashed 
lines show the accepted kinematic parameters PA and i. The error-bars on the V ro t(r) curve correspond to the mean 
non-circular velocities (see the text for details). 

Fig. 6. The same as figEJfor NGC 23. 

Fig. 7. The same as figEJfor NGC 615. 

Fig. 8. The same as figEJfor NGC 972. 

Fig. 9. The same as figEJfor NGC 1084. 

Fig. 10. The same as figEJfor NGC 1134. 

Fig. 11. The same as figGjJfor NGC 2964. 

Fig. 12. The same as figElfor NGC 3583. There is a ghost-image in the bottom-left corner of the monochromatic and 
continuum images. 



Fig. 13. The same as figEJfor NGC 3893. 
Fig. 14. The same as figCOfor NGC 4100. 
Fig. 15. The same as figEJfor NGC 4136. 
Fig. 16. The same as figCOfor NGC 4414. 
Fig. 17. The same as figEJfor NGC 4814. 
Fig. 18. The same as figEJfor NGC 5371. 



Fig. 19. The same as figOOfor NGC 6643. 
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